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Abstract. We consider a Fabry-Perot resonance (a transmission band edge 
resonance) in periodic layered structures involving birefringent layers. In our 
previous publication [2] we have shown that the presence of birefringent layers with 
misaligned in-plane anisotropy can dramatically enhance the performance of the 
photonic-crystal resonator. It allows to reduce its size by an order of magnitude 
without compromising on its performance. The key characteristic of the enhanced 
slow-wave resonator is that the Bloch dispersion relation uj (k) of the periodic 
structure displays a degenerate photonic band edge, in the vicinity of which the 
dispersion curve can be approximated as Auj ~ (Afc)*, rather than Au) ~ (Afc)^. 
Such a situation can be realized in specially arranged stacks of misaligned 
anisotropic layers. On the down side, the presence of birefringent layers results 
in the slow wave resonance being coupled only with one (elliptic) polarization 
component of the incident wave, while the other polarization component is 
reflected back to space. In this paper we show how a small modification 
of the periodic layered array can solve the above fundamental problem and 
provide a perfect impedance match regardless of the incident wave polarization, 
while preserving the giant slow-wave resonance characteristic of a degenerate 
photonic band edge. Both features are of critical importance for many practical 
applications, such as the enhancement of various light-matter interactions, light 
amplification and lasing, optical and microwave filters, antennas, etc. 
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1. Introduction 

Wave propagation in spatially periodic media, such as photonic crystals, can be 
qualitatively different from any uniform substance. The differences are particularly 
pronounced when the wavelength is comparable to the primitive translation L of the 
periodic structure O IH [5l [6l [71 [H [9] . The effects of strong spatial dispersion culminate 
at stationary points ujs = to (ks) of the Bloch dispersion relation where the group 
velocity u = doj/dk of a traveling Bloch wave vanishes 

duj 

— = 0, at k = ks, uj = uJs = to {ks) . (1) 

One reason for this is that vanishing group velocity always implies a dramatic increase 
in density of modes at the respective frequency. In addition, vanishing group velocity 
also implies certain qualitative changes in the eigenmode structure, which can be 
accompanied by some spectacular effects in wave propagation. A particular example of 
the kind is the frozen mode regime associated with a dramatic amplitude enhancement 
of the wave transmitted to the periodic medium [TUl El [121 [HI Ull [E]- In this 
paper, we focus on a different slow-wave effect, namely, on a Fabry-Perot resonance 
in bounded photonic crystals. This slow wave phenomenon, illustrated in Figs. [H 
and m is also referred to as the transmission band edge resonance. There are some 
similarities between the frozen mode regime and the slow-wave resonance in plane- 
parallel photonic crystals. Both effects are associated with vanishing group velocity at 
stationary point of the Bloch dispersion relation. As a consequence, both effects 
are strongly dependent on specific type of spectral singularity ([1]). A fundamental 
difference though is that the frozen mode regime is not a resonance phenomenon in 
a sense that it is not particularly sensitive to the shape and size of the photonic 
crystal. For instance, the frozen mode regime can occur even in a semi-infinite 
periodic structure, where the incident plane wave is converted to a frozen mode slowly 
propagating through the periodic medium until it is absorbed [lOl [Til [12 [ISl [HI [H] • 
By contrast, in the case of a slow wave resonance, the entire bounded periodic structure 
acts as a resonator, resulting in a strong sensitivity of the resonance behavior to the 
size and shape of the photonic crystal. 

It is also important to distinguish between two qualitatively different classes of 
photonic-crystal resonators. The first class comprises resonance cavities where the 
role of periodic dielectric structure reduces to electromagnetic (EM) field confinement 
by reflecting it back to the cavity interior. The resonance frequency (or frequencies) 
of such photonic cavities usually lies in a frequency gap (a stop-band) of the photonic 
crystal. The periodic dielectric array here plays the role of a distributed Bragg 
reflector. The number of resonance modes depends on the cavity size. It can be 
a single mode localized on an isolated defect inside the photonic crystal [T71 [Hj ■ Or 
the cavity can support multiple resonances, if its size significantly exceeds the light 
wavelength. More detailed information on photonic crystal cavities can be found in 
numerous papers and monographs on optics and photonics (see, for example, [19j and 
references therein). In this paper, we will not further discuss this subject. 

The second class of photonic-crystal resonators comprises the, so-called, slow- wave 
resonators. They are qualitatively different from the band-gap cavities. In slow-wave 
photonic- crystal resonators, the reflectors may not be needed at all, as shown in the 
example in Fig. [H The role of the periodic structure here is to support slow EM 
waves. The resonance frequencies lie in the transmission bands of the photonic crystal 
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- not in band gaps. Although, the resonance frequencies can be very close to a band 
edge (dJ, as shown in Fig. [21 A typical example of slow wave resonance in a photonic 
crystal is presented by the transmission band edge resonance, illustrated in Figs. [T] 
and [21 In certain cases, slow- wave resonators can provide significant advantages over 
cavity resonators. They are used for the enhancement of light-matter interactions, 
such as nonlinear and nonreciprocal effects, optical activity, light amplification and 
lasing, etc. They can also be used in optical and microwave filters, delay lines, 
as well as for the enhancement of antenna gain and directionality. More detailed 
information can be found in an extensive literature on the subject (see, for example, 
[20l[6l[7l[i[2Tl[22l[23l[22[25],and references therein) . 

In this paper we describe a slow-wave photonic-crystal resonator with drastically 
reduced dimensions and enhanced performance, compared to that of a common Fabry- 
Perot resonator based on a periodic stack of non-birefringent layers. The idea is to 
employ periodic structures supporting the dispersion relations different from those 
allowed in periodic arrays of non-birefringent layers. Indeed, periodic arrays involving 
birefringent layers can display stationary points (|T]) different from a regular photonic 
band edge in Fig. [2ta). Some examples are shown in Fig. [3l Slow waves associated 
with such stationary points can produce giant transmission band-edge resonances, 
much more powerful compared to those achievable in common layered structures. The 
first step in this direction was made in 2] , where it was shown that the transmission 
band-edge resonance in the vicinity of a degenerate photonic band edge (DBE) in 
Fig. [3jb) produces much better results, compared to a regular photonic band edge 
(RBE) of Figs. [2ja) and[3ja). Specifically, at the frequency of DBE related giant 
slow-wave resonance, the electromagnetic energy density inside the photonic-crystal 
can be estimated as 

(Wdbe) oc WiN\ (2) 

where W/ is the energy density of the incident wave and N is the total number of 
unit cells in the periodic stack. By comparison, the average EM energy density at a 
regular transmission band-edge resonance in Fig. [2|is 

(Wrbe) (X WiN^. (3) 

The estimations ([2]) and ([3]) imply that the Q-factor of a DBE based slow-wave 
resonator can be by factor N"^ higher compared to that of a RBE related Fabry- 
Perot resonator of the same size. And this a huge difference! A detailed comparative 
analysis of the giant DBE related slow wave resonance versus the regular transmission 
band edge resonance can be found in [2]. 

On the down side, periodic structures with birefringent layers have a fundamental 
problem - their reflectance and transmittance are essentially dependent on the 
incident wave polarization. This dependence is particularly strong near the edges 
of transmission bands, where the slow-wave resonances occur. In particular, a DBE 
related giant transmission resonance described in [2] is coupled only with one (elliptic) 
polarization component of the incident wave, while the other polarization component 
is reflected back to space by the photonic crystal boundary [2]. In other words, at the 
resonance frequency, a periodic stack involving birefringent layers acts as a polarizer, 
reflecting back to space roughly half of the incident wave energy. This behavior is 
illustrated in Figs. [Hand [5] Similar problem exists in all different modifications of the 
frozen mode regime considered in [TUl [HI [HI [T31 [HI [H] . For many applications, such a 
polarization selectivity may not be acceptable. In this paper we offer a solution to the 
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above problem. We show how to utihze aU the incident wave energy, while preserving 
the extraordinary performance of the DBE based photonic-crystal resonator. The 
idea is to modify the periodic layered array so that instead of a degenerate band edge, 
the respective dispersion curve develops a split photonic band edge (SBE) shown in 
Fig. [3jb). Under certain conditions specified below, the photonic resonator with 
a SBE will display a giant transmission band edge resonance, similar to that of a 
DBE. But, in addition, the SBE resonator couples with the incident wave regardless 
of its polarization and, therefore, utilizes all the incident EM radiation - not just one 
polarization component. The latter feature is of critical importance for a variety of 
practical applications. 

Similar approach can be applied not only to a photonic-crystal cavity resonance, 
but also to all different modifications of the frozen mode regime described in 

[iniiiiiiisiiiiiiii]. 



2. Dispersion relation and energy flux in periodic arrays of birefringent 
layers 

In this section we introduce some basic definitions and notations of electrodynamics of 
stratified media composed of birefringent layers. A detailed and consistent description 
of the subject, along with numerous references, can be found in [2l [lOl [TTl [15] , where 
similar notations and terminology are used. For simplicity, we restrict ourselves to 
the case of a plane monochromatic wave normally incident on a layered structure, as 
shown in Fig. [T] The results can be easily generalized to the case of oblique incidence, 
as it was done, for example, in 115| . 



2.1. Transverse waves in stratified media 



Time-harmonic electromagnetic field inside and outside the layered medium can be 
described by the vector-column 



E^{z) 

Ey{z) 
Hy{z) 



(4) 



where E [z) and H (z) are time-harmonic electric and magnetic fields. The z direction 
is normal to the layers. The values of ^ at any two different locations z and z' are 
related by the transfer matrix T {z, z') defined by 



^(z) ^T{z,z')-^ {z') . 



(5) 



The elements of the transfer matrix are expressed in terms of material tensors and 
other physical characteristics of the stratified medium. 

Let ^E*/, "ifii, and be the incident, reflected, and transmitted waves, 
respectively, as shown in Fig. [T]. To the left of the stack (at z < 0), the electromagnetic 
field is a superposition of the incident and reflected waves. To the right of the stack 
(at z > D), there is only the transmitted wave. The field inside the periodic medium 
is denoted as ^t- Since all four transverse field components in ([4]) are continuous 
functions of z, we have the following boundary conditions at z = and z = D in Fig. 

m 



(0) + (0) = *T (0) , 



(6) 
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*P (D) = {D) . (7) 
2.2. Time-harmonic field inside periodic layered structure 

At any given frequency to, the time- harmonic field (z) inside the periodic stratified 
medium can be represented as a superposition of Bloch eigenmodes, each of which 
satisfies the following relation 

*fc(^ + i) ^e^'^'^^fc (8) 

or, equivalently, 

VP, (z) = e^'^^Vfe (z) , (z + i) = V^fe (z) . (9) 

The Bloch wave number k is defined up to a multiple of 2t:/L. The correspondence 
between cu and k is referred to as the Bloch dispersion relation. Real k correspond to 
propagating (travehng) Bloch modes. Propagating modes belong to different spectral 
branches to (k) separated by frequency gaps. In reciprocal and/or centrosymmetric 
periodic structures, the Bloch dispersion relation is always symmetric with respect to 
the points k = and k = n/L oi the Brillouin zone 

uj{ka + k) =Lu{ko-k), (10) 

where 

fco = 0, tt/L. (11) 

In periodic structures composed of non-birefringent layers, every Bloch wave is 
doubly degenerate with respect to polarization. Typical k — oj diagram for such a case 
is shown in Fig. [2l^a). If, on the other hand, some of the layers display an in-plane 
anisotropy or gyrotropy, the polarization degeneracy can be lifted. The respective 
k — uj diagrams are shown in Fig. [3) 

Finally, if some of the layers are magnetically polarized and, in addition, the 
periodic array is non-centrosymmetric, the dispersion relation can also develop spectral 
asymmetry 

UJ {k) ^ UJ {~k) . (12) 

This effect can be significant if the magnetic layers display an appreciable 
nonreciprocal circular birefringence (magnetic Faraday rotation). Examples of 
periodic layered structures with asymmetric dispersion relation (fT2| can be found 
in [TBI [To]. Further in this paper we only discuss non-magnetic structures. 

The speed of a traveling wave in a periodic medium is determined by the group 
velocity u = dw/dk. Normally, every spectral branch w (fc) develops stationary points 
([T} where the group velocity of the corresponding propagating mode vanishes. Usually, 
such points are located at the center and at the boundary of the Brillouin zone 

fc, = fco = 0, tt/L. (13) 

This is always the case in periodic layered structures composed of non-birefringent 
layers, where all stationary points coincide with photonic band edges, as shown in 
Fig. [IJa). If, on the other hand, some of the layers in a unit cell are birefringcnt, 
then in addition to (fT3|) . some dispersion curves can also develop a reciprocal pair of 
stationary points corresponding to a general value of the Bloch wave number fc, as 
shown in Fig. [3jb). The respective portion of the k — u) diagram can be described as 
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a split band edge (SBE). The dispersion relation can develop a DBE or a SBE only if 
the periodic layered array has birefringent layers with misaligned in-plane anisotropy 
[21 [TS]. Example of such a layered structure is shown in Fig. [HI 

Unlike propagating modes, evanescent Bloch modes are characterized by complex 
wavenumbers k ^ k' + ik" . Under normal circumstances, evanescent modes decay 
exponentially with the distance from the periodic structure boundaries. In such cases, 
the evanescent contribution to \1/t can be significant only in close proximity of the 
surface or some other defects of the periodic structure. The situation can change 
dramatically in the vicinity of a stationary point H]) of the dispersion relation. At 
first sight, stationary points ([T]) relate only to propagating Bloch modes. But in fact, in 
the vicinity of every stationary point frequency LOg, the imaginary part k" of the Bloch 
wavenumber of at least one of the evanescent modes also vanishes. As a consequence, 
the respective evanescent mode decays very slowly, and its role may extend far beyond 
the photonic crystal boundary. In addition, in some special cases, the electromagnetic 
field distribution in the coexisting propagating and/or evanescent eigenmodes becomes 
very similar, as u approaches LUg. This can result in a spectacular effect of coherent 
interference, such as the frozen mode regime [101 HH El HI]- What exactly happens 
in the vicinity of a particular stationary point essentially depends on its character and 
appears to be very different for different types of singularity (fl]). 

Using the transfer matrix ([5]), the Bloch relation ([5]) can be recast as 

TL*fc = e^'^^k (14) 

where is the transfer matrix of a unit cell 

Tl^T{L,0), *fe = *fe(0). 

At any given frequency, there are four Bloch eigenmodes, propagating and/or 
evanescent, with different polarizations and wave numbers 

(z) , *fe2 (z) , ^-fca (z) , *fc4 (z) ■ (15) 

Depending on the frequency w, the full set (jlSp of Bloch eigenmodes may include only 
propagating modes, only evanescent modes, or both. In any event, the respective set 
[fci, k2,k^, fc4] of four Bloch wave numbers satisfies the relation 

[fci, fc2, fcg, = [K.k*, k;,kl] , (16) 

which is a direct consequence of J unitarity of the transfer matrix T (see [15j and 
references therein). 

Taking into account (jl6p . one can distinguish the following three possibilities. 

(i) All four Bloch modes in are propagating 

ki = kl, fc2 = ^2, k^ ^ k^, ki — kX (17) 

In this case, two of the propagating modes have positive group velocities, they 
are referred to as forward waves. The other two Bloch waves have negative group 
velocities, they are referred to as backward waves. 

(ii) All four Bloch modes in are evanescent 



ki 7^ fci, fc2 7^ fcj, fcs 7^ kl, ki ^ kX 



(18) 
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This is the case when the frequency to falls into a photonic band gap. According 
to one can assume that 

ki = fcj, fc3 = kl- (19) 

Two of the evanescent modes have k" > 0; they are referred to as forward 
evanescent modes. The other two evanescent modes have k" < 0; they are referred 
to as backward evanescent modes. 

(iii) Two of the Bloch modes in ([15]) are propagating modes, while the other two are 
evanescent. 

fci = kl, k2 = k;, fca ^ k;, ki ^ kl (20) 
Again, according to (jl6p . one can assume that 

kl = kl, k2 = fcjj = K- (21) 

In all cases, propagating modes with m > and evanescent modes with k" > 
are referred to as forward waves. The propagating modes with u < and evanescent 
modes with k" < are referred to as backward waves. 



2. 3. EM energy flux in layered media 

The real- valued energy flux (the Poynting vector) associated with time-harmonic field 
© is 



S 



ReE{z) X ReH{z) 



(22) 



The expression (|22]) can also be recast in the following compact form [l4] 



where 



J 



5=-(vI/,JvI'), 



1 
0-10 
0-100 
10 



(23) 



(24) 



2.3.1. EM energy flux in layered media In the case of a lossless stratified medium, 
the Poynting vector 5* in (|23)) is independent of the coordinate z. 

In a periodic layered medium, the column vector \E' (z) in (f23|) can be represented 
as a linear combination of four Bloch components ((T5|) 

■^T {z) = *fei {z) + *fe2 [z) + *fe3 {z) + ^-M {z),0<z<D. (25) 
Substituting ^ into ^ yields 

1 4 

At this point we can apply the orthogonality relation 

(^-^J^-jO^O, iffc, ^fc*, (27) 
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the proof of which can be found in [T3], P. 327. The expression (|26p together with 
(|27p lead to the following conclusions regarding the energy flux of a time-harmonic 
electromagnetic field in lossless periodic layered media. 

1) The contribution of each propagating Bloch mode to the total energy flux is 
independent of the presence or absence of other Bloch eigenmodes with the same 
frequency 

^ = E^^ = ^E(*"'^*»)' (28) 

1=1 i=l 

where the summation runs over all propagating eigenmodes. The number of 
propagating modes can be 4, 2, or 0, depending on which of the cases pT|) . (PO)) . 
or (fT5|) we are dealing with. 

2) The contribution of evanescent eigenmodes to the energy flux depends on how 
many of them exist at this particular frequency lo. 

(a) In the case pO| of two evanescent modes ^3 and ^"4 we have 

5 = Re (*3, , where = fcg, (29) 

which implies that only a pair of evanescent modes with conjugate wave 
numbers can contribute to the energy flux. The respective contribution 
([^ is independent of the presence of propagating modes ^'i and '^2- In 
accordance with Eq. (|29p . a single evanescent mode, either "if 3 or ^4, does 
not produce energy flux on its own. 

(b) In the case ([T8| of four evanescent modes we have 

S = Re J*2) + Re (*3, J*4) , where fca = fc*, k^ = k*, (30) 

which implies that either of the two pairs of evanescent modes with conjugate 
wave numbers contribute to the energy flux independently of each other. 

In the general case of the time-harmonic EM field (z) being a superposition 
(|25p of several Bloch modes, propagating and evanescent, the contribution of each 
propagating mode ^'fc to the total energy flux can be expressed in terms of its group 
velocity and amplitude 

Sk = WkUk, (31) 

where 

W,^{\^,\)' = {\^,\f (32) 

In the particular case of a single propagating mode, the quantity Wk in ()3ip is equal 
to the energy density averaged over a unit cell L. 

By contrast, a single evanescent mode never transfers the energy. Only the 
combination (|29p or pop of forward and backward evanescent modes can contribute 
to the energy flux. 



2.3.2. Transmission/reflection coefficients Let Si, Sr, St, and Sp be the energy 
fluxes of the respective waves in Fig. [TJ The transmission/reflection coefficients t and 
r are defined as 

Sp Sp .„„^ 
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In the case of losses medium, the Poynting vector S is independent of z 



S — Si + Sr = St = Sp. 



In such a case 



t = l-r = S/Si. 



(34) 



Otherwise, 



t + r + a 



1, 



where a > is the absorption coefficient of the stack. 
2.4- Single mode excitation regime 

To explain the nature of strong polarization dependence of the transmission band edge 
resonance in periodic stacks of anisotropic layers we introduce some basic concepts 
related to semi-infinite periodic layered structures. Then, we will see how those 
concepts apply to finite periodic stacks. 

Consider a plane monochromatic wave incident on a semi-infinite periodic 
structure with the same unite cell L as that of the finite periodic stack in Fig. [T] 
In other words, let us assume that the right-hand boundary of the photonic crystal in 
Fig. [T]is absent. In such a case, the EM field (z) inside the periodic medium is 
a superposition of two forward Bloch waves, each of which can be either propagating 
with w > 0, or evanescent with k" > 0. The backward Bloch modes (propagating with 
u < and evanescent with k" < 0) do not contribute to (z) in the case of a semi- 
infinite periodic structure. According to 15J, the incident wave polarization can always 
be chosen so that only one of the two forward Bloch waves is excited. Such a situation 
is referred to as a single mode excitation regime. The respective two polarizations 
of the incident wave are referred to as the single mode excitation polarizations. In 
periodic structures involving birefringent layers, a single mode excitation polarization 
is, generally, elliptic. Note, that if all the layers are non-birefringent, then, due 
to polarization degeneracy, the very concept of a single mode excitation regime at 
normal propagation makes no sense. In the case of oblique incidence, though, a single 
mode excitation regime exists even in periodic structures composed of non-birefringent 
layers; the respective incident wave polarizations are linear and referred to as TE and 



As we already mentioned, the dispersion relation of a periodic structure involving 
birefringent layers usually does not display polarization degeneracy. In particular, the 
periodic stack in Fig. [6] can display polarization degeneracy at normal propagation 
only if the misalignment angle (|87|) between the adjacent anisotropic layers Ai and A2 
is 7r/2. In all numerical examples in this paper, the misalignment angle (/) is different 
from and 7r/2, because otherwise the existence of a DBE and a SBE would be ruled 
out by symmetry 15J. Therefore, none of the spectral branches in Fig. [3]is degenerate 
with respect to polarization. 

2.4.1- The vicinity of a RBE Consider the vicinity of a RBE gi or 52 in Fig. [3l[a). If 
the incident wave frequency lies in a transmission band close enough to the respective 
RBE, the transmitted wave {z) inside the semi-infinite periodic layered structure 
is a superposition of two forward Bloch waves 



TM. 



*T {z) = -^pr {z) + -^ev {z) , 



z > 



(35) 
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one of which is propagating and the other is evanescent. There are two single mode 
excitation regimes at such a frequency. The first one produces only propagating mode 
^pr (z) at z > 0, while the other produces only evanescent mode 5'e„ (z) at z > 0. 
The respective transmission coefficients satisfy the following relations 

if (z) = ^pr (z) , then t = tpr > 0, (36) 

if (z) = *et, (z) , then t = tev = 0. (37) 

Eq. ((37)) reflects the fact that according to (129]) . a single evanescent mode does not 
contribute to the energy flux. 

The relations ([M]) and ([57)1 imply that only one of the two polarization 
components of the incident wave (the one producing the forward propagating mode 
^pr) can be transmitted to the semi-infinite photonic crystal. The other polarization 
component of the incident wave (the one producing only evanescent mode 5'et,) will 
be totally reflected back to space. 

Now, what happens if the periodic structure has a finite thickness? In such 
a case, the relation ([35]) does not apply. Instead, the EM field (z) inside the 
finite periodic stack is a superposition of all four Bloch waves ([T5|) . including both 
forward and backward modes. If the number N of unit cells in the periodic stack is 
significant, then at certain frequencies, the pair of propagating modes (one forward 
and one backward ) form a standing wave (|40p constituting the transmission band 
edge resonance. An important point, though, is that such a resonance is excited only 
by one (elliptic) polarization component of the incident wave - the same component 
that produces a single propagating mode inside the semi-infinite structure. If, on the 
other hand, the incident wave polarization is such that it produces a single evanescent 
mode inside the semi-infinite structure, then it will be reflected back to space from the 
finite periodic stack as well. In practice, the relation (l37|) is valid even for the finite 
periodic structures composed of as few as several unit cells L. 

2.4-2. The vicinity of a DBE Consider now the vicinity of a DBE d in Fig. [Hb). If 
the incident wave frequency lies in the upper transmission band close to the DBE 
d, the transmitted wave ^'t (z) inside the semi-infinite periodic layered structure 
still is a superposition (j35p of one propagating and one evanescent Bloch waves. 
But remarkably, in close proximity of the DBE, the column vectors ^'p^ and in 
(j35p become nearly parallel to each other, while their amplitude diverges [15]. More 
specifically, at the photonic crystal boundary at z = 

^pr (0) « -fe^ (0) OC |CJ - CJrfP^/* , US UJ ^ UJd, (38) 

provided that the incident wave has some general polarization and fixed amplitude. 
The destructive interference ([55]) ensures that the boundary condition ^ can be 
satisfied, while both Bloch contributions to \1/t (z) diverge. As the distance z from 
the slab boundary increases, the evanescent component (z) dies out 

■^ev (z) - "fev (0) exp i-zk") , (39) 

while the propagating component "ilpr (z) remains constant and very large. Eventually, 
as the distance z further increases, the transmitted wave (z) reaches its large 
saturation value corresponding to its propagating component ^!pr (z). This constitutes 
the DBE related frozen mode regime in a semi- infinite photonic crystal |15] . 
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Generally, the incident wave polarization can always be chosen so that only one 
of the two forward Bloch waves in (j35p is excited in the semi-infinite photonic crystal 
(a single mode excitation regime). What happens to the frozen mode regime in such 
a case? Obviously, the destructive interference ([55)1 cannot occur in a single mode 
excitation regime. As a consequence, the incident wave is reflected back to space 
without producing the frozen mode, regardless of whether the single mode excitation 
regime produces a single propagating or a single evanescent Bloch mode |15j . The 
latter circumstance is related to the fact that the incident wave polarizations producing 
a single propagating mode and a single evanescent mode become indistinguishable as 

The fact that all four Bloch modes ([15]) are nearly parallel to each other in the 
vicinity of a DBE is also very important for understanding the DBE related giant 
transmission band edge resonance in a finite photonic crystal. Indeed, as we pointed 
out earlier, in the case of RBE related resonance, the incident wave producing a 
single propagating mode is coupled with the resonance mode, while the incident wave 
producing a single evanescent mode is not. By contrast, in the case of DBE related 
resonance, both polarizations are nearly indistinguishable. As a result, neither of the 
two single mode regime polarizations is coupled with the resonance mode. Instead, 
the best result is achieved if the incident wave polarization is orthogonal to that of 
single mode regime polarization. A numerical illustration of this point is given in Fig. 

m 

In summary, the DBE related transmission band edge resonance is not excited if 
the incident wave polarization corresponds to the single mode excitation regime. This 
implies that if the incident wave polarization is randomly chosen, then on average, 
at least half of the incident wave energy is reflected back to space from the photonic 
crystal boundary, without contributing to the resonance. The above problem can 
be overcome in slightly modified periodic structures displaying a SBE, rather than a 
DBE. 

3. Simplified geometrical description of transmission band-edge resonance 

Steady-state Fabry-Perot resonance in a pane-parallel photonic crystal is commonly 
described as a standing wave composed of a pair of reciprocal propagating Bloch waves 
(|40p with equal and opposite Bloch wave numbers and vanishing group velocities. The 
nodes of the standing wave coincide with the photonic crystal boundaries at z = 
and z — D. This simple representation works very well in the vicinity of a regular 
photonic band edge (RBE). But, in the cases of DBE or SBE related giant transmission 
band edge resonance, the above simple picture should be significantly modified. We 
start this section with a brief description of the standing wave concept as applied to a 
regular transmission band edge resonance. Then we proceed to a comparative analysis 
of giant transmission band edge resonances in the vicinity of a DBE and SBE. 

3.1. Slow-wave resonance as a standing wave 

Consider a finite plane-parallel periodic structure composed of N unit cells L. If 
the number N is significant, the electromagnetic properties of the periodic stack can 
be effectively described in terms of the Bloch eigenmodes ([9]) of the corresponding 
infinite periodic layered medium. For instance, a periodic stack having as few as 
several periods L can display almost total reflectivity at photonic band gap frequencies, 
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which resuhs from the absence of propagating Bloch modes at the respective frequency 
range. A typical dependence of the finite stack transmittance on the incident wave 
frequency is shown in Fig. [21Ib). The sharp peaks in transmission bands correspond 
to transmission band edge resonances. The resonance frequencies are located close 
to a photonic band edge, where the group velocity is very low. At each resonance, 
the electromagnetic field (2) inside the periodic stack is close to a standing wave 
composed of a pair of reciprocal Bloch waves with large and nearly equal amplitudes 
and low group velocity 

(z) = ^-fc (z) + *-fe (z) , (40) 



where 



k^k^, < z < D = NL. 



The photonic crystal boundaries at z = and z — D coincide with the standing wave 
nodes, where the forward and backward Bloch components interfere destructively to 
meet the boundary conditions ^ and Q. The latter circumstance determines the 
wave numbers of the Bloch eigenmodes contributing to the resonance field inside the 
periodic structure 

TT 

{kjn - ko) ^ ±—m, m = l,2,..., (41) 

The wave number ko is given in (fTTI) . The integer m denotes the resonance peaks 
in order of their distance from the respective photonic band edge. The respective 
resonance frequencies are expressed in terms of the dispersion relation 

w„=w(fc™). (42) 

According to (|4ip . the proximity of the resonances to the photonic band edge is 
determined by the number N of unit cells in the periodic stack. The expressions 
(|4T|) and only apply if 

TV > m. (43) 

Let us introduce the following dimcnsionless notations for the small deviation of 
the wave number and the frequency from the respective stationary point 

fi^\k-ko\L, i'{k) ^[u;{k)~Lj{ko)]L/c. (44) 

According to (|4ip and (|^^ . the resonance values of k and ly are 

TT 

Km ~ -^TO < 1, Vyn~V (k„) . (45) 

The most powerful resonance is usually the one closest to the respective photonic band 
edge 

Ki ~ < 1, vi Kv{ni). (46) 
3. 2. Energy density and energy flux at resonance frequency 

Assume that the amplitude of the incident wave in Fig. [l]is unity. Due to the boundary 
condition ([6]), we have at the left-hand boundary of the photonic crystal in Fig. [1] 



(0) = -^k (0) + -^-k (0) oc 1. 



(47) 
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At resonance frequency, the transmission coefficient (I34p is close to unity 

i„ = S/Si oc 1. (48) 

Therefore, the relation similar to ([T7|) also applies at the right-hand boundary of the 
periodic stack 

^TiD) ^■^k{D) + ^-k{D) o^l. (49) 

Let us see under what conditions the finite resonance transmission (|48p is 
compatible with the vanishing group velocity in the vicinity of stationary point ([1]), 
where the transmission band edge resonance occurs. According to pi[) . the energy 
flux associated with the time-harmonic field (z) in (|40)) is 

S = Sk + S-k'^ Wku (k) + W-kU i-k) = [Wk ~ W-k) u [k) , (50) 

where 

u{>,)^\u{k)\ = \u{-k)\, (51) 

and 

In the vicinity of stationary point ((T|), the magnitude u (k) of the group velocity in 
(|50p is vanishingly small. The fact that the resonance energy flux remains of the order 
of unity implies that the amplitude of the Bloch components in ((50)) increases so that 

Wk - W-k oc as M ^ 0, (52) 

or, equivalently, 

(l^-fel") - cx u-\ &su-^Q. (53) 

Thus, finite transmittance at m — s- requires the divergence (|53p of the amplitude of 
the Bloch components. 

The requirements (|47)) and (|49l) together with (|53p imply that although the 
amplitudes of the forward and backward Bloch components both diverge as m ^ 0, 
they must interfere destructively near the photonic crystal boundaries 

*fc (0) « (0) cx 

and 

In order to reconcile the boundary condition (|47p with the requirement (j53p of a finite 
energy flux, we have to impose the following requirement on the amplitudes of the two 

Bloch components 

Wk - W^k oc ^/Wk oc u-^, as m ^ 0. (54) 

The relation (|54p was derived under the assumption that the time-harmonic filed 
{z) inside the periodic medium is a superposition ([40]) of one forward and one 
backward propagating Bloch waves with equal and opposite wave numbers and group 
velocities. This assumption is warranted in the vicinity of a regular photonic band 
edge of a periodic stack with birefringent layers, such as those shown in Fig. [S^a). 
In a periodic stack of non-birefringent layers, at any given frequency there are two 
pairs of reciprocal Bloch waves with equal wave numbers but different polarizations. 
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Still, because of the polarization degeneracy, the above consideration can be directly 
applied to the latter well. 

By contrast, in the DBE and SEE related transmission band edge resonances, the 
contribution of all four Bloch components to the resonance field 'i'r (z) is absolutely 
essential. In such cases, the treatment of the transmission band edge resonance as a 
simple standing wave (|50p is not adequate. 

3.3. Regular photonic band edge (RBE) 



Let us apply the relation (|54p to the classical case of Fabry-Perot resonance in 
the vicinity of a regular photonic band edge (RBE). Using the notations ((44|l . the 
dispersion relation in the vicinity of a RBE can be approximated as 



(55) 



The group velocity magnitude is 



u ~ OKKi y'2v/a, Q < u/a. (56) 

According to (j45p . the resonance values of the wave number and the frequency are, 
respectively 

7rm a /Trmy 



2 



The group velocity magnitude at resonance is 



Substitution of Um to ([54]) yields 



am 

The energy density distribution W {z) at RBE related resonances is typical of a 
standing wave 

W{z)ozWi(—\ sin^f— z), m = l,2,.... (58) 
\m J \ NL I 

In the case of the most powerful first transmission resonance ( the one with 
m = 1), the two Bloch components in ([^0]) interfere constructively in the middle of 
the periodic structure. The respective field intensity is 

cx WiN"^. (59) 
3.Jf-. Degenerate photonic band edge (DBE) 

The dispersion relation in the vicinity of a DBE d in Fig. [3^b) can be approximated 
as 

,.{k)^^^\ (60) 

The magnitude u of the group velocity of each of the reciprocal pair of propagating 
modes is 

w « « V8 {bv^) < biy. (61) 
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Unlike the case of RBE related resonance, in the vicinity of a DBE, the EM field inside 
the periodic medium cannot be approximated by a simple standing wave (gOl) • Instead, 
we have to take into account the evanescent mode contribution, which appears to be 
comparable to that of the propagating modes [Ij. As the result, the simple geometrical 
conditions (|4T|) and (l45l) give only rough order of magnitude estimate of the resonance 
values of the wave number and frequency. More importantly, the maximum field 
intensity inside the photonic crystal at the frequency of the transmission band-edge 
resonance is now 

Wk oc WiN'^, (62) 

which is by factor N'^ larger than that of a regular transmission resonance given in 
dSS]). 

The detailed analysis of the DBE related giant transmission resonance was carried 
out in [2]. Here, we would like to highlight the down side of this effect, which is the 
polarization selectivity illustrated in Figs. S] and [5l As we already stated before, 
the DBE related giant transmission resonance is coupled only with one (elliptic) 
polarization component of the incident wave, while the other polarization component 
is reflected back to space by the photonic crystal boundary. The reflected polarization 
component corresponds to the single mode excitation regime, explained at the end of 
Section 2. A solution to this fundamental problem is to slightly modify the periodic 
structure so that instead of a DBE, the respective spectral branch displays a SBE. 
Example of a SBE is shown in Fig. [Hb). 

4. Transmission band-edge resonance in the vicinity of a split photonic 
band edge (SBE) 

Simplified description of SBE related resonance 

In the previous section we described the resonance conditions for a single pair (|40p 
of reciprocal Bloch waves. In periodic stacks involving birefringent layers, a slow- 
wave resonance produced by a single pair of reciprocal Bloch waves can be excited by 
only one polarization component of the incident plane wave. The other polarization 
component is not coupled with the respective resonance mode. This implies that if 
the incident wave polarization is random, at least half of the incident wave energy will 
not enter the periodic medium and will be reflected back to space by the photonic 
crystal boundary. 

Consider now the vicinity of a split photonic band edge (SBE) on the k — cu 
diagram in Fig. ^h). The physical characteristics of the periodic structure are chose 
so that the SBE in Fig. [Hb) is close to a DBE. The frequency range 

oja < U! < ujb, (63) 

covers a narrow portion of the transmission band which includes the SBE. At any 
given frequency from (|63p . there are two pairs of reciprocal Bloch waves with very 
low group velocity and different polarizations; each pair being capable of producing 
its own slow-wave cavity resonance with the resonance conditions similar to (I4ip or 
(|45p . Our focus is on the possibility of the two resonances occurring at the same 
frequency. Such a situation will be referred to as the double resonance. It turns out 
that the double transmission band edge resonance is as powerful as the giant resonance 
associated with a DBE. But, in addition, a SBE related resonance utilizes all the 
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energy of the incident wave regardless of its polarization. By contrast, a DBE based 
giant transmission resonance is coupled only with one polarization component of the 
incident wave; the rest of the incident wave energy being reflected back to space. This 
important difference is obvious if we compare the DBE related transmission dispersion 
shown Fig. |4]and the SBE related transmission dispersion shown in Fig. [T] 

To start with, let us consider the dispersion curve with a SBE in more detail. If 
the split between the twin band edges bi and 62 in Figs. Eljb) is small, the dispersion 
relation in the vicinity of SBE can be approximated as 

i^in)^ln^+^-^\ (64) 



where 
and 



a/b < (65) 



\a/b\ < 1. (66) 
The inequality (p5|) is the condition for SBE. Indeed, in the opposite case of 

a/b > 0, (67) 

the dispersion curve (IM)) would develop a RBE at k = 0, as shown in Fig. El^a). While 
in the case 

a/b = 0, (68) 

the dispersion curve would have a DBE at k = 0. The additional inequality 
is the condition for the proximity of the SBE to a DBE. This proximity allows us to 
use the expansion (|64p in the frequency range spanning both twin edges of the SBE. 
More importantly, the condition ([55]) is essential for the phenomenon of the giant 
transmission resonance in the vicinity of SBE. 

There are three stationary points associated with a SBE. The first one is trivial 

Ka - 0, l^a= 0. (69) 

It is located either at the center of the Brillouin zone, or at its boundary. The other 
two stationary points correspond to the actual SBE 



±Kb = ±yj-alb, = -a^/Ah. (70) 

Taking into account ([701) . condition for the proximity of the SBE to a DBE 
can be recast as 

Kb « 1. (71) 

The condition (|7ip implies that the points 61 and 62 on the dispersion curve are close 
to each other. 

In what follows, we assume for simplicity that 

6 < < a. (72) 

In this case, the SBE in question corresponds to the upper edge of the transmission 
band, as shown in Fig. [IJb). The alternative case of 



a < < 5 



(73) 
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corresponds to the SBE being the lower edge of the respective transmission band. 
There is no quaUtative difference between the two cases. 
At any given frequency v within the range 

I^a<l^ < Vb, (74) 

there are two pairs of reciprocal Bloch waves. Each pair comprises one forward and 
one backward propagating modes with equal and opposite wave numbers and group 
velocities 



±Kh^/l - Jl - Va<V<Vb. (75) 



± Kes = ±K6 Jl + - -JJ-, l'<Vb, (76) 



The pair of wave numbers ([75|) corresponds to the concave portion of the dispersion 
curve (|64|) , while the pair of wave numbers (|76|) corresponds to the convex portion of 
the dispersion curve. Obviously, 

Km < Kex, at Va <V < Vb- (77) 

At any given frequency, the EM field inside the periodic stack is 

{z) = {Z) + ^-e. [z] , (78) 

where 

{z) = ^',„. (z) + (z) , (79) 

and 

{z) = ^^^^ {z) + (z) . (80) 

To distinguish between 5'i„ (z) and "^ex {z)^ we will refer to the respective quantities 
as the external and internal. The inequality (|77p justifies these terms. 



4-2. Conditions for the double SBE resonance 

Within the frequency range ([74|) . either pair of the reciprocal Bloch waves ([79]) and 
(|80p can develop a transmission resonance. Of particular interest here is the case 
where the two resonances occurs at the same or almost the same frequency. This 
situation is illustrated in Fig. [Hljc), as well as in Fig. [3 and [8] 

Let us start with the internal resonance created by the reciprocal pair ([79|l of Bloch 
waves corresponding to the concave section of the dispersion curve. It is possible that 
the frequency range (Tfl]) contains only a single cavity resonance - the one with m = 1. 
Such a case is determined by either the proximity of the SBE to a DBF, or by the 
right choice of the number N of the unit cells in the stack. According to ((70)) and 
(|4ip. the condition for a single internal resonance is 

TT 

Ki < Kfc < 2ki, = ki = — . (81) 
The respective resonance frequency Vin is determined by (|42| and ([64]) 



a 2 ^4 



(82) 
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where ki = n/N. 

Consider now the external resonance created by the reciprocal pair ([50]) of Bloch 
waves corresponding to the convex section of the dispersion curve. Let us impose the 
condition 

I'ex = l^in = Vr (83) 

that the external and internal resonances occur at the same frequency (f82|) . This 
condition leads to the following equality 

2 b 4 0-2 ^4 

Simple analysis shows that it is only possible if at = i^r we have 

_2tt 

The relation 

^ex 2^^271 

together with and yield 

/5" /5"7r 
The frequency of the double transmission resonance is 




The group velocities of the two reciprocal pairs of Bloch waves at the resonance 
frequency Vr are 

and ^ 

We. = ±36^3 = ±35 (^) . 

By comparison, in the case of a DBE related giant transmission band edge 
resonance, we would have the following estimation for the resonance frequency vi 
and the respective group velocity of the two Bloch components 

b b 

These estimations are similar those related to the double transmission resonance 
associated with SBE. In either case, the average resonance energy density is estimated 
by ([3]), which justifies the term 171071^ transmission resonance. 

Let us remark that the entire consideration of this subsection was based on the 
assumption that each pair ([75]) and (|80p of the reciprocal Bloch waves is responsible 
for its own individual transmission resonance, described as a standing wave (j79p or 
(|80p. respectively. While the double resonance at is described as the situation 
where the frequencies Vin and v^x of those individual resonances merely coincide. In 
fact, these two transmission resonances can be treated as independent only if the 
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respective resonance frequencies a-nd I'ex are well separated. As soon as and 
Vex are close to each other, the EM field '^t {z) in ([75]) becomes a superposition of all 
four propagating Bloch modes. The latter situation persists even if the incident wave 
polarization correspond to the, so-called, single mode excitation regime, defined in the 
previous section. In other words, the single mode excitation regime produces almost 
pure internal or external resonance, provided that their frequencies are well separated. 
Otherwise, the EM field '^t (z) is a superposition of all four Bloch eigenmodes with 
two essentially different sets of wave numbers. In such a case, the resonance field 
4't (^) cannot be viewed as a standing wave (201) • The physical reason for such a 
strong hybridization is that due to the condition (|7T|) . the REE in question is very 
close to a DBE. On the other hand, in the vicinity of a DBE, all four vector-columns 
(z) in ([25|) become nearly parallel to each other [SlII^. The latter circumstance 
excludes the possibility of exciting only one of the two pairs of the reciprocal Bloch 
modes ([75]) or ([50]) in the situation (|83p . where the resonance conditions (|¥T|l are 
in place for both of them simultaneously. The strong hybridization implies that the 
estimation (I84p of the double resonance frequency is rather crude. Still, it provides a 
very useful guidance on the conditions for SBE related giant transmission resonance. 
The bottom line is that the SBE related giant transmission resonance is as powerful as 
that related to a DBE. But, in addition, the SBE related resonance provides a perfect 
coupling with the incident wave regardless of its polarization. 

5. Conclusion 

In summary, we would like to stress that the remarkable features of the DBE and 
SBE related giant transmission resonances can be derived from such fundamental 
characteristics of the periodic composite medium as its electromagnetic dispersion 
relation. Specific details of the periodic array, such as physical characteristics of the 
constitutive components, or structural geometry, are only important as long as the 
symmetry of the periodic array is compatible with the existence of the required spectral 
singularities. 
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6. Appendix 



The simplest periodic layered structure supporting a DBE or a SBE at normal 
propagation is shown in Fig. \E\ A unit cell L contains one isotropic B layer and 
two misaligned anisotropic layers Ai and A2 with inplane anisotropy. The isotropic 
layers have the thickness B and the dielectric permittivity 



£b 



SB 
es 









(85) 



The dielectric permittivity tensors £a in each anisotropic A layer has the form 



+ (5 COS 2(y9 (5sin2(/9 
5 sin 2ip £a^5 cos 2<y9 

£3 



(86) 



where the parameter 5 characterizes the magnitude of inplane anisotropy, and the 
angle determines the orientation of the anisotropy axes in the x — y plane. All the 
A layers have the same thickness A and the same magnitude 5 of inplane anisotropy. 
The only difference between the adjacent anisotropic layers Ai and A2 in [6] is their 
orientation ip. An important characteristic of the periodic structure in Fig. [6] is the 
misalignment angle 

between the layers Ai and ^2- This angle determines the symmetry of the periodic 
array and, eventually, the kind oik — u diagram it can display. 

In all numerical simulations related to the periodic layered structure in Fig. [6] we 
use the following values of the material parameters in (|85p. ([86|l and ([87|) 



£b^ 16.0, EA = 4.7797, 5 = 3.4572, = 7r/6. 



At normal propagation, the numerical value of £3 in (j86p is irrelevant. The relative 
thickness of the A and B layers, can be different in different examples. 

In all plots of the field distribution inside periodic media ai < z < D we, in 
fact, plotted the following physical quantity 



(89) 



\^{z)\') = (E{z)-E* {z) + H{z)-H* (z) 



which is the squared field amplitude averaged over a local unit cell. The actual function 
I'l' (z)l^, as well as the electromagnetic energy density distribution W (z), are strongly 
oscillating functions of the coordinate z, with the period of oscillations coinciding 
with the unit cell length L. Given the relation oc l^t^ (z)]^, the quantity ((89| can 
also be qualitatively interpreted as the smoothed energy density distribution, with the 
correction coefficient of the order of unity. 

In all plots, the distance z, the wave number k, and the frequency w are expressed 
in units of L, and cL~^, respectively. 
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Figure 1. Scattering problem of a plane wave normally incident on a periodic 
stack of dielectric layers. The indices R, and P denote the incident, reflected 
and transmitted waves, respectively. The field inside the periodic medium is 
In the case of a slow wave resonance, the incident wave frequency lies in a 
transmission band of the periodic structure, close to a band edge, as illustrated 
in Fig. [2] 
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Figure 2. (a) A fragment of a typical Bloch k — u) diagram of a periodic array 
composed of non-birefringent layers; gi and g2 are the edges of the lowest photonic 
band gap. Each spectral branch is doubly degenerate with respect to the wave 
polarization, (b) Transmission dispersion t (oj) of the respective finite periodic 
stack; the sharp peaks near the edges of the transmission bands are associated 
with slow-wave Fabry-Perot resonances, also known as transmission band edge 
resonances. The location I I42I I of the resonance peaks depends on the number A'^ 
of unit cells L in the periodic stack. 
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Figure 3. Fragments of the k — u! diagrams of the periodic layered structure in 
Fig. [6] for two different values of the ratio A/B of the layer thicknesses. Regular, 
degenerate , and split photonic band edges are denoted by symbols g, d, and b, 
respectively. The Bloch wave number k and the frequency lu are expressed in 
units of 1/L and c/L. 
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Figure 4. Transmission dispersion of the periodic stack of 18 unit cells at 
frequency range including a DBE at = oj^j. The respective k — lu diagram 
is shown in Fig. |3jb). In the cases (a) and (b), the incident wave in linearly 
polarized. In the cases (c) and (d), the incident wave polarization is adjusted so 
that at any given frequency it corresponds to a single mode excitation regime: in 
the case (c) it is a single propagating mode, while in the case (d) it is a single 
evanescent mode. Obviously, in the latter case the incident wave is reflected back 
to space. 
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Figure 5. Smoothed energy density distribution W (z) at frequency of the first 
(closest to the DBE) giant transmission band edge resonance in Fig. |4]for four 
different polarizations of the incident wave. In a single mode excitation regime of 
Figs, (c) and (d), the transmission resonance is suppressed. Particularly so in the 
case (d), where the EM field inside the periodic medium corresponds to a single 
evanescent mode. 
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Figure 6. Periodic layered structure with a unit cell L containing two misaligned 
anisotropic A layers, and one isotropic B layer. The respective dielectric 
permittivity tensors are given in I I85I I , II86I I , and 11881 . This is the simplest layered 
array supporting the Bloch dispersion relation with a DBE and/or a SBE, as 
shown in Fig. [3Tb). 
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Figure 7. Manifestation of SBE related double resonance in the transmission 
dispersion t (ui) of periodic stack with N = 18. The respective k — to diagram 
is shown in Fig. [3b). Observe that at the resonance frequency, the stack 
transmittance is close to unity regardless of the incident wave polarization. By 
contrast, in the case of DBE-related giant transmission resonance in Fig. |4] the 
impedance matching is polarization dependent. In the cases (c) and (d), the 
incident wave polarization is adjusted so that at any given frequency it would 
correspond to the respective single mode excitation regime in the semi-infinite 
layered structure. 



Slow wave resonance in periodic stacks of anisotropic layers 



29 




Figure 8. Smoothed energy density distribution at frequency of the SBE related 
giant transmission resonance in Fig. [7] for four different polarizations of the 
incident wave. The cases (c) and (d) relate to a single mode excitation regime. 
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Figure 9. Transmission dispersion of periodic stacks composed of different 
number A'^ of unit cells. The frequency range shown includes SBE on the k — u) 
diagram in Fig. [3b). The red and the blue curves correspond to two different 
polarizations of incident wave - in either case, the incident wave polarization is 
adjusted so that at any given frequency it corresponds to a single mode excitation 
regime. In the case (b) oi N = 18, the two resonance frequencies nearly coincide, 
creating condition for double transmission resonance with perfect impedance 
matching. 



